www.gusucode.com > 循环自相关函数工具箱源码程序 > matlab代做 修改 程序循环自相关函数工具箱/cyclostationary_toolbox/ex6.m
clear; clc; % t=0:1/600:1023/600; % A=1;B=1;fb=7;fa=100; % x=A*(1+B*cos(2*pi*fb.*t)).*cos(2*pi*fa.*t); % figure(1) % plot(t,x) % alpha=2*pi/300; load('E:\169.mat'); %%读取数据 % speed=1796; %motor rotational speed(rpm) sig1=X169_DE_time'; ln=length(sig1); n=1000; t=0:1/12000:(n-1)/12000; figure(1) plot(t,sig1(1:n)) xlabel('\fontsize{9}\fontname{标楷体}时间/s');ylabel('\fontsize{9}\fontname{标楷体}幅值/um'); max_tau=0; alpha=2*pi/12000; R=cyclic_autocorrelation(sig1(1:n),alpha,max_tau); y=1000*R.*conj(R); % ff=600*(0:1023)/512; ff=0:1:6000-1; figure(2) % plot(ff(1:150)/(2*pi),y(1:150)); plot(ff(1:3000),y(1:3000)); xlabel('\fontsize{9}\fontname{标楷体}循环频率 /Hz');ylabel('\fontsize{9}{R_{x}}^{\alpha}(0) /um^2');